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ABSTRACT 


A  five  degree-of-freedom  simulation  has  been  constructed  for  each  of  four 
candidate  tactical  weapons  for  use  in  a  program  to  develop/evaluate  new 
techniques  for  emitter  homing.  The  vehicle  types  include  a  powered  missile,  a 
projectile,  a  minidrone,  and  a  glide  bomb.  The  roll  degree-of-freedom  was 
eliminated  from  the  simulation  on  the  basis  that  each  vehicle  is  stabilized  In 
roll  with  low  roll  rates.  The  report  provides  detailed  descriptions  of  the 
simulations  Including  the  aerodynamic  and  autopilot  models,  along  with  the 
numerical  procedures  used.  Typical  responses  of  the  closed  loop  homing 
performance  are  given  for  each  vehicle. 
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1.0  INTRODUCTION 


This  report  describes  a  five  degree-of-freedom  computer  model  of  four 
different  vehicle  types.  These  models  are  intended  for  use  in  closed  loop 
simulations  to  develop/evaluate  new  techniques  for  emitter  homing  [1,]  and/or 
other  autonomous  seekers. 

The  four  vehicles  modeled  are:  1)  a  powered  missile,  2)  a  cannon- launched 
guided  projectile,  3)  a  mini-drone,  and  4)  a  glide  bomb.  The  powered  missile 
is  based  upon  dynamic  and  aerodynamic  characteristics  that  would  be  typical  of 
the  Maverick.  The  projectile  model  is  based  on  parameters  associated  with 
the  Copperhead  cannon-launched  guided  projectile  [2,  6,  7,  8].  Preliminary 
data  from  the  XBQM-106  was  used  to  develop  the  mini-drone  model  [9],  Finally, 
the  GBU-20  (GBU-15  PWW)  [10,  11,  12]  was  the  basis  for  the  glide  bomb.  It  is 
emphasized  that  while  the  simulations  are  considered  realistic,  they  should 
not  be  considered  as  exact  simulations  for  these  specific  vehicles. 

In  all  of  these  models,  the  roll  degree-of-freedom  has  been  eliminated  by 
assuming  that  the  vehicle  is  roll  stabilized  with  relatively  low  roll  rates. 
In  the  case  of  the  mini-drone,  "side  force"  surfaces  were  Incorporated  to 
enhance  the  yaw  maneuvering  capability  without  the  additional  complexity  of 
bank-to-turn  maneuvers. 

Section  2.0  provides  a  detailed  description  of  the  coordinate  system,  the 
equations  of  motion,  the  aerodynamic  forces,  the  line-of-sight  (LOS) 
sensor/processor  and  the  guidance/autopilot  model  used  for  each  vehicle 
type.  Section  3.0  outlines  the  difference  equations  used  for  the  simulations 
that  were  developed  from  the  data  in  Section  2.0.  Finally,  Section  4.0 
presents  typical  output  from  the  simulation  as  mechanized. 


2.0  MODEL  DESCRIPTION 


2.1  General 

Figure  2.1  Illustrates  the  overall  geometry  of  the  vehicle  and  its  local, 
or  body,  coordinate  system  with  respect  to  the  fixed  space  coordinate 
system.  Upper  case  letters  (X,  Y,  Z)  are  used  to  denote  the  fixed  coordinate 
system  and  lower  case  letters  used  for  the  body  fixed  coordinate  system. 
Figure  2.2  is  a  companion  figure  used  to  further  illustrate  and  define  the 
geometrical  parameters  and  their  relationship  to  one  another. 

Since  the  vehicle  is  roll  stabilized  only  two  angles,  the  Euler  angles  0 
and  are  needed  to  define  the  attitude  of  the  vehicle  body  axes  with  respect 
to  the  space  axis  system.  It  is  further  noted  that  the  body  axes  are  rotated 
with  respect  to  the  fixed  axes  so  that  for  8  and  »  0,  the  vehicle  x-axis  is 
-X,  and  the  vehicle  z-axis  is  -Z. 

The  Q  and  the  R  represent  body,  angular  rotation  rates  (pitch  and  yaw) 
about  the  y  and  z  axes,  respectively.  It  is  further  noted  that 

0  -  /  Qdt  (2.1a) 

♦  -/S-jdt  (2.1b) 

The  velocity  of  the  vehicle  center  of  mass  (the  origin  of  the  body  axis) 
is  defined  by  the  three  velocity  components  u,  v,  w,  along  the  x,  y,  z  axes, 
respectively. 

•  •  • 

From  Fig.  2.1,  the  velocity  components,  X^,  Y  ,  Z^  of  the  vehicle  center 
of  mass  are  given  by 

Xv  *  -u  co8  0  cosi|;  +  v  sin<|»  -  w  8in0  cos i|/  (2.2a) 

Yv  -  u  cos0  sin4>  +  v  co8^  +  w  sin0  siniji  (2.2b) 

Zy  -  u  8in0  -  w  co80  (2.2c) 
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Fiq.  2.1.  Geometrical  relationships 


An  overall  block  diagram  of  the  simulation  is  shown  in  Fig.  2.3.  As  may 
be  seen,  there  are  four  main  blocks  or  models.  The  sensor  feedback  between 
the  airframe  and  the  autopilot/guidance  blocks  is  shown  in  a  general  sense, 
that  is,  each  model  does  not  necessarily  Include  all  of  the  sensors  as 
shown.  The  different  models  will  utilize  various  combinations  as  will  be  seen 
below.  The  feedback  paths  without  sensors  are  natural  or  inherent  feedback 
due  to  the  governing  equations.  For  example,  there  is  an  altitude  feedback 
path  to  the  airframe  model.  Its  function  is  to  supply  the  information  for  the 
altitude  dependent  aerodynamic  parameters.  In  the  case  of  Q  and  R  (pitch  and 
yaw  body  angular  rates)  both  the  natural  and  the  sensor  feedback  paths  are 
used  from  the  airframe  to  the  LOS  sensor/processor  model.  The  reason  for  this 
is  that  the  proportional  navigation  model  used  requires  that  the  angle  between 
the  body  axis  and  the  LOS  must  be  subtracted  from  the  sensor  measurement  since 
a  strapdown  sensor  has  been  assumed.  This  will  be  discussed  further  when  the 
model  details  are  considered. 

The  kinematic  block  of  the  simulation  is  given  by  equations  2.1  and  2.2 
along  with  the  following  additional  integration 

X  -  /  X  dt 

v  4  v 

Y  -  /  Y  dt 

v  ‘  v 

Z  -  /  Z  dt 

V  1  V 

The  blocks  representing  airframe  and  the  guidance  autopilot  models  will  be 
considered  in  some  detail  in  the  following  subsections.  An  elementary 
consideration  of  the  LOS  sensor/processor  will  be  given  here.  A  more  detailed 
description  is  given  in  reference  2.  The  description  given  below  will  not 
discuss  the  details  of  how  the  measurement  is  made  with  the  antenna  sensor 
which  interacts  in  a  complex  way  with  the  target  and  terrain. 


llOMTHl 


Fig.  2.3.  Overall  block  diagram  of  simulation  model. 
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2.2  Airframe  Equations  of  Motion 

The  equations  of  motion  of  a  rigid  body  with  five  degress  of  freedom  (the 
roll  degree  of  freedom  is  not  included)  are  well  known  (Refs.  3,  4)  and  are  as 
follows: 


Q  -  My 

R  -  M 
z 

u  *  A  -  g  sln8  -  —  Qw+  —  Rv  +  T  /m 
x  r  r  x 

o  o 

v  «  A  -  —  Ru  +  T  /m 
y  r  y 

o 

w  ■  A  +  g  cos0  +  —  Qu  +  T  /m 
z  r  z 

o 

where  rQ  -  180/x ,  the  dot  over  the  variable  indicates  differential  with 
respect  to  time,  and 


Q  »  body  pitch  rate  (in  deg  per  sec) 

R  ■  body  yaw  rate  (in  deg  per  sec) 

u  *  body  axial  velocity  (along  x  direction  -  ft/sec) 
v  -  body  lateral  velocity  (along  y  direction  -  ft/sec) 
w  -  body  vertical  velocity  (along  z  direction  -  ft/sec) 

My  ■  normalized  aerodynamic  moment  about  Y  axis  or  pitching  moment 

Mz  “  normalized  yaw  aerodynamic  moment 

A^j  -  normalized  aerodynamic  axial  force 

Ay  »  normalized  aerodynamic  sideforce  (in  yaw  plane) 

Az  *  normalized  aerodynamic  normal  force  (in  pitch  plane) 

Ty ,  Ty ,  Tz  *  x,  y,  z  components  of  the  thrust  in  lbs. 

(NOTE:  that  if  the  thrust  is  misaligned  with  respect  to  the  center  of  mass, 

additional  moment  contributions  will  be  induced.)  The  normalized  moments  and 
forces  are  given  as  follows: 


7 


M  -  C  - =: 

y  m  I 


q  S  c  r 

_ _ y  ° 


yy 


m  -  c 

z  n 


q  S  c  r 
z  o 


zz 


A  -  C  (Si) 
x  x  * 


A  -  C  (Si) 
y  y  vm  ' 


A  -  C  (Si) 
*  z  vm  ' 
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where  q  -  ^  pv  ■  dynamic  pressure  (lbs/ft  ) 

2  2  2 

v  ■  u  +  v  +  w  ,  total  vehicle  velocity  (ft/sec) 


p  ■  air  density  at  specified  altitude  (or  z)  (slugs  per  cubic  foot) 


S  «  reference  area  (ft  ) 


m  ■  vehicle  mass  ■  W/g  (slugs) 
w  «■  vehicle  weight  (lbs) 

2 

g  ■  gravity  acceleration  -  32.3  ft/sec 

c ,  c  ■  reference  lengths  associated  with  the  y  and  z  directions, 
y  respectively  (ft) 

2 

1,1  »  pitch  and  yaw  moments  of  inertia,  respectively  (slug  ft  ) 

yy  zz 

The  C  ,  C  ,  C  ,  C  ,  C  are  aerodynamic  coefficients  which  are  functions,  in 
tn  n  x  y  z 

general,  of  the  angles  of  attack  in  pitch  and  yaw  (a,  0  respectively)  as  well 
as  Mach  no  (M)  and  the  pitch  and  yaw  control  surface  deflections  (6g,  6^), 
respectively.  The  Mach  no.  is  defined  as 
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-  i 


M  -  V/a 

o 

where  aQ  is  the  speed  of  sound  and  is  a  function  of  altitude. 

The  values  of  the  aerodynamic  coefficients  are  usually  determined  by  wind 
tunnel  measurements,  although  analytical  estimates  are  sometimes  used. 
Depending  upon  the  vehicle  characteristics  and  velocity  range,  the 
coefficients,  range  from  simple  to  complex  requiring  extensive  two  or  three 
dimensional  tables. 

It  is  further  to  be  noted  that  some  of  the  inertial  parameters  can  vary 
with  time  because  of  rocket  motor  bum  off.  For  convenience,  however,  the 
pertinent  values  of  the  steady  state  parameters  of  all  of  the  vehicles  are 
given  in  Table  2.1. 


2.3  Aerodynamic  Forces 

Up  to  this  point,  the  formulation  is  general  and  applies  to  airframe 
model  of  all  of  the  vehicles.  The  aerodynamic  coefficients  must  be 
specialized  for  each  of  the  vehicles  which  is  done  in  appendices  as  follows: 


Appendix  1 
Appendix  2 
Appendix  3 
Appendix  4 


Missile 
Projectile 
Minidrone 
Glide  Bomb 


The  data  in  these  appendices  provides  the  aerodynamic  coefficients  as 
functions  of  M,  a,  8,  6g,  <5R.  In  addition,  there  are  time  varying  inertial 
parameters  for  the  Missile  to  account  for  the  change  of  mass,  moment  inertia, 
and  center  of  mass  as  the  rocket  motor  burns.  Note  that  this  occurs  in  the 
first  four  seconds. 

The  variation  of  the  speed  of  sound  and  air  density  is  given  by  the 
following  analytical  expressions: 

p  -  Po  (1  -  1.11  Z  +  0.36  Z2) 


9 


I 


TABLE  2.1 

GEOMETRIC  AMD  INERTIAL  PARAMETERS 
(STEADY  STATE  VALUES) 


PARAMETER 

MISSILE 

PROJECTILE 

MINIDRONE 

GLIDE  BOMB 

S  (FT2) 

0.786 

0.196 

18.18 

16.6 

W  (LBS) 

398 

135 

135 

3000 

1^  (SLUG  FT2) 

53.5 

6.13 

25.5 

646.2 

Izz  (SLUG  FT2) 

53.5 

6.13 

32.4 

716.6 

S 

1 

0.5 

1.65 

1.54 

c„ 

1 

0.5 

1.65 

11.33 

p  -  0.002378 
o 

Z  -  Z/ 40000 

a  -  1117  -  (4.22  x  10_3)Z  0  <  Z  <  36000 

o 

a  -  965  36000  <  Z 

o 

where  p  is  in  slugs/ft J  (  *  lb-sec  Vft^)  and  aQ  is  in  ft/sec.  Other  formulas 
or  tables  could  be  used,  but  these  provide  a  reasonable  variation  of  the 
parameters  needed  for  trajectory  calculations. 

2.4  LOS  Sensor/Processor 

As  noted  in  Fig.  2.3,  the  inputs  to  the  sensor/processor  are  the  position 
in  space  of  the  vehicle  and  target,  along  with  the  actual  and  measured  body 
rates  in  pitch  and  yaw.  Basically,  the  task  of  the  sensor/processor  is  to 
measure  the  line-of-sight  between  the  vehicle  and  the  target  and  from  this 
measurement  provide  a  command  to  the  guidance  computer.  Since  the  sensor  is 
fixed  to  the  body,  the  body  angle  must  be  subtracted  from  the  actual  sensor 
measurement  to  obtain  the  LOS.  This  is  done  with  the  body  rate  sensors. 
Since  the  rate  sensors,  or  gyros,  are  not  perfect,  the  errors  associated  with 
them  must  be  included  in  the  simulation  model.  In  addition,  the  sensor  Itself 
will  introduce  inaccuracies  which  also  must  be  included.  The  actual 
simulation  of  the  sensor  measurement  is  extensive  including  reflections, 
multipath,  terrain  irregularities,  and  signal  interference  will  not  be 
considered  here. 

The  true  LOS  in  elevation  and  azimuth  are  given  as 
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where  the  subscripts  Indicate  vehicle  and  target.  To  these  LOS  angles  must  be 
added  the  errors  as  noted  above.  The  gyro  errors  in  elevation  and  azimuth,  6£ 
and  H>e>  respecitvely ,  are 


e  -  e 


cN 


0  n  +  0  „ 
eB  eR 


♦e  "  *eN  +  *EB  +  *eR 

where  the  subscripts  N,  B,  and  R  indicate  noise,  bias,  and  receiver, 
respectively.  The  noise  and  bias  components  from  the  rate  gyro  can  be 
expressed  as 


eN  " 

i «« « 

eB  ’ 

I  0.,  1* 

R 

1  «s 

f  JSLdt 

eN 

J  cos0 

R 

f  _£B_ 

J  CO80 

eB 

where  the  cosO  is  needed  to  transform  the  body  yaw  angular  rate  to  the  azimuth 
plane.  This  bias  value  of  pitch  and  yaw  rate  (QeB,  ReB)  varies  from  sensor  to 
sensor  and  from  one  period  of  time  to  another.  The  values  are  usually  given 
as  lo  level  or  a  maximum  not  to  exceed  level  in  degrees/sec  or  degrees/hr. 
The  noise  component  of  the  error  is,  in  general,  correlated  and  can  be 
expressed  as  [14] 


+  °G 


n 


2_ 


n 


where  tg  is  the  correlation  time  of  the  noise 
Oq  is  the  rms  level  of  the  noise 


n  is  unit  intensity  white  noise 


The  spectrum  for  this  model  of  the  noise  is 


<Kw)  “ 


2  tc°g 

1  +  (T  G(0); 


where  u>  is  the  circular  frequency  in  radians  per  second.  For  simulation 
purposes,  the  unit  intensity  white  noise  level  is  generated  from  a  random 
sequence  (or  random  number)  generator.  This  will  be  discussed  further  in 
Section  3.0. 

Thus,  the  measured  LOS  can  be  expressed  as 


AZ«  •  \  *  8e»  ♦  8eB  +  8e* 


YM  Y  eN  eB  TeR 


The  measured  values  are  then  filtered,  differentiated  and  then  filtered 
again.  For  example,  the  measured  elevation  LOS  is  processed  as  follows: 


FiHer:  dfc  <XZMF>  ”  "  Tf  (XZMF5  +  Tf  XZM 


Differentiate 
Numerically:  X 


XZMF  -  XZMF  , 
n  n-1 


ZMF  (t  -  t  ,) 
n  n-1 ' 


Filter:  dt  <XZMFF)  “  "  Tf  (XZMFF^  +  Tf  ^ZMF5 


•  1  • 

—  (X. 

ZMFFF'  Tf  v  ! 
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The  frequency  at  which  this  is  done  can  be  substantially  slower  than  the 

general  Integration  of  the  airframe  equations  of  motion.  The  fact  that 

several  Integrations  of  the  system  can  occur  over  the  sample  time  (tn  -  tn_j) 

is,  in  effect,  a  sample  and  hold.  The  value  X_.__  provides  the  proportional 

ZMFF 

navigation  signal  to  the  guidance  computer  in  pitch  or  elevation.  The  value 
determines,  in  part,  the  switching  to  terminal  guidance  for  the  glide 

ZMr  P  F 

bomb. 


The  measured  pitch  “gimbal”  angle  is  defined  as 


9 CM  XZM  ”  9 


As  for  the  line-of-sight,  the  measured  sensor  gimbal  angle  is  filtered  before 
being  used  as 

dt  (0gmP  "  “  9GM1  +  9gm 

d /  n  \  _  _  ^  a  4,  i o 

dt  '  GM2;  t_  °GM2  T  r  GM1 

B  D 

The  filtered  values  of  the  gimbal  angle  are  used  for  both  mode  switching  and 
for  trajectory  shaping  as  indicated  below. 


2.5  Guidance/Autopilot 

As  for  the  aerodynamic  forces,  the  guidance/autopilot  model  must  be 
specialized  for  each  vehicle  type.  The  function  of  this  model,  as  suggested 
by  the  subsection  title,  it  twofold.  The  autopilot  function  is  to  insure  that 
the  vehicle  is  stable  and  controllable.  Each  of  the  four  vehicle  types 
considered  requires  some  supplementary  stability  augmentation  to  insure 
adequate  response  over  the  complete  flight  regime.  The  guidance  requirement 
is  to  convert  the  two  commands  from  the  LOS  sensor/processor  into  appropriate 
positions  of  the  control  surfaces  to  guide  the  vehicle  to  the  target.  This 
part  of  the  model  is  also  used  to  represent  any  control  surface  dynamics  and 
limits  that  are  appropriate. 
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Figures  2.4,  2.5,  2.6,  and  2.7  are  models  of  the  guidance/autopilot 
models  for  the  missile,  projectile,  minidrone,  and  glide  bomb,  respectively. 
As  may  be  seen,  there  is  a  fairly  wide  variation  from  model  to  model. 
Generally,  the  autopilots  are  more  complex  than  the  models  given  here.  The 
simplifications  are  based  primarily  on  eliminating  the  higher  frequency 
components.  For  example,  the  rate  gyros  would  probably  have  frequency 
components  greater  than  20  Hz.  Thus,  it,  along  with  any  shaping  and  filter 
currents,  has  been  replaced  by  pure  gains.  The  basic  underlying  assumption 
here  is  that  adequate  provision  has  been  made  to  insure  the  stability  of  the 
autopilot.  When  this  has  been  done  properly,  these  higher  frequency 
components  are  not  needed  for  the  guidance  studies  for  which  these  models  are 
designed.  The  parameters  not  given  on  the  diagrams  are  given  in  Table  2.2. 
The  mode  switch  sequences  for  the  glide  bomb  are  given  in  Table  2.3.  The 
programmed  yaw  gain  schedule  given  in  Table  2.3  reduces  the  yaw  gain  during 
the  activation  of  the  pitch  acceleration  autopilot  in  the  terminal  phase  of 
flight.  The  reason  for  this  is  to  suppress  the  undesirable  yaw/pitch  coupling 
at  high  angles  of  attack. 

As  may  be  seen,  the  inputs  for  all  of  the  models  are  the  filtered  line- 
of-sight  rate  in  both  pitch  and  yaw  (azimuth  and  elevation)  as  well  as  the 
effective  gimbal  angle  which  is  used  to  shape  the  trajectory.  The  output  from 
the  autopilot  is  the  control  surface  deflections  in  pitch  and  yaw.  The 
elevation  and  azimuth  navigation  gains,  GZN  and  Gyjg.  have  been  made  "noise 
adaptive"  to  suppress  the  deleterious  effects  of  interference  "noise"  produced 
in  a  multiple  emitter  environment.  This  "noise"  generally  produces  excessive 
maneuvers  early  in  the  trajectory  before  "locking  on"  to  a  specific  target. 
These  excessive  maneuvers  cause,  in  turn,  a  substantial  reduction  in  velocity, 
and  hence  range,  thereby  decreasing  the  accuracy  markedly.  The  approach  here 
is  to  reduce  the  gain  under  these  severe  circumstances.  The  form  of  the  gain 
is 


G 


where 


i-N  - 

V  -  l  (X  -  X)2 

i“N-n 

o 
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Fig.  2,4.  Missile  guldance/autopllot  model 
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TABLE  2.2 

VALUES  FOR  CONSTANTS  USED  IN  SIMULATION 


UNITS 

MISSILE 

PROJECTILE 

MINIDRONE 

GLIDE  BOMB 

coz 

g/DEG/SEC 

3.5 

10.47 

2  DEG/DEG/SEC 

3 

COY 

g/DEG/SEC 

3.5 

9.13 

2  DEG/DEG/SEC 

2 

*A 

DEG/g 

2.2 

0.8 

— 

— 

*Q 

DEG/DEG/SEC 

0.2 

0.15 

0.15 

1 

*R 

DEG/DEC/SEC 

0.2 

0.15 

0.15 

1 

T1 

SEC 

0.03 

— 

— 

— 

t2 

SEC 

0.10 

— 

— 

— 

t3 

SEC 

0.05 

— 

0.1 

— 

*G 

g/DEG 

0.0175 

0.35 

0 

— 

l/rf 

Hz 

0.33/0.25* 

0.33/0.25 

0.25/0.17 

0.33/0.25 

1/tA,B 

Hz 

0.33/0.50 

0.33/0.5 

0.25/0.25 

0.33/0.50 

Kzl 

DEG/g 

— 

— 

— 

1.25 

*YL 

DEG/G 

— 

— 

— 

2.50 

kzb 

DEG/DEG/SEC 

— 

— 

— 

0.25 

kyb 

DEC/DEG/SEC 

— 

— 

— 

0.25 

*2 

DEG/SEC/g 

\ 

— 

— 

5 

DEG/SEC/g 

— 

— 

5 

kzh 

DEG/SEC/DEG 

— 

— 

— 

1.0 

^YH 

DEG/SEC/DEG 

— 

— 

— 

1.0 

ao 

DEG 

— 

— 

— 

5.5/-12.0 

N 

— 

— 

— 

— 

2.00 

“ref 

DEG 

— 

— 

— 

4 

8  REF 

DEG 

20 

28.6 

20 

20 

8c 

DEG 

-17.2 

-11.50 

** 

-20 

°BIAS 

8 

3.5 

4.25 

— 

1.0 

nzl 

8 

— 

— 

— 

2.5 

"yl 

8 

1.5 

PITCH/YAV 

‘Activated  after  4  second*  of  flight 


ill 


TABLE  2.3 


MODE  SWITCH  SEQUENCE  FOR  GLIDE  BOMB 


1.  Flight  Initiated  with  a-Hold  node 

2.  Glmbal  hold  occurs  when  the  following  conditions  are  met 


2.2  Flight  time  >  3  seconds 

2.3  9  <  9.7° 

3.  Acceleration  auutopllot  (terminal)  occurs  when  the  following  conditions 
are  met 

3'1  *ZMfFF  <  -°‘344  °/SEC 

3.2  Flight  time  >  3  seconds 

3.3  8  <  -5.73° 

3.4  Global  hold  oust  have  been  activated 

4.  ip  command  not  used  for  the  trajectories  considered 
3.  Taw  gain  schedule 


5.1 

gyn  “  2,0 

0  *  t  S  C glmbal  hold 

5.2 

gyn  “  0,75 

Cglmbal  hold  ^  C  *  terminal  +  10 

5,3 

gyn  -  l'n 

ttcrmlaal  <  t 
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N  -  current  tine  interval 


nQ  -  number  of  sample  in  the  window 


A^  ■  ith  value  of  the  LOS  rate 

—  ,  i-N 

A  *  —  2  A. ,  or  the  average  value  of  the  LOS  rate  in  the  window 

1  no  i-N-n  1 

o 

Vo  -  reference  variance  of  the  LOS  rate 


Gq  »  nominal  noise  free  gain 


For  periods  of  very  high  noise  levels,  this  law  is  quite  effective  in  reducing 
the  loss  in  energy  caused  by  the  noise  induced  maneuver  drag. 

It  is  to  be  noted  that  for  some  vehicles  (notably  the  glide  bomb  and  the 
projectile)  that  the  use  of  the  noise  adjusted  gain  is  not  effective  for  the 
yaw  channel  and  is,  therefore,  not  used.  The  primary  reason  for  this  is  that 
for  sluggish  vehicles  (low  maVieuver  capability),  there  is  not  sufficient 
excursion  for  target  acquisition. 

This  variable  gain  can  be  implemented  by  means  of  recursion  filters  as 


vi  “  x2  "  xl 

i  i 


where  V  is  the  variance  at  the  current  Instant,  and 


x,  -  (l-A)  x«  +  AA, 
zi  zi-l  1 


x.  *  (l-A)  x.  +  AA. 
*i  *1-1  1 


A  -  exp(-T/ty) 


where  tv  is  the  filter  time  constant 
T  is  the  simulation  Interval 
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3.0  SIMULATION  EQUATIONS  AND  PROCEDURE 


The  sequence  of  Che  simulation  proceeds  In  the  following  steps: 

1.  Enter  Initial  conditions 

2.  Update  all  variables 

3.  Guidance/autopilot  commands 

4.  Aerodynamic  forces 

5.  Equations  of  motion 

6.  Kinematics 

7 .  LOS  processor 

8.  Return  to  Step  2 

Each  one  of  these  steps  will  be  outlined  and  the  required  difference  equations 
will  be  given. 

3.1  Initial  Conditions 

The  Initial  conditions  for  the  simulation  are  as  follows: 
a.  The  position  of  the  vehicle  and  target: 

*V  *T 

o  o 


b.  The  velocity,  heading  and  angle-of -attack  of  the  vehicle 


The  velocity  and  angles-of-attack  are  used  to  determine  the  body 
velocity  components  as 

u  ■  V  >/l  -  tan^a  -  tan^J 
oo  o  o 

v  -  u  tang 
o  o  o 
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w  •  u  tana 
o  o  o 

From  these  values  the  Initial  LOS  and  global  angles  can  be  computed 
using  the  equations  provided  in  Section  2.4. 

c.  The  control  surface  positions 

6Eo*  6Ro 

d.  All  of  the  initial  values  of  the  filters  are  used  in  the 

simulation.  Those  values  not  specified  are  taken  to  be  zero.  In  the 
case  of  the  gyro  noise  which  will  be  discussed  further  in  Seclon  3.6, 
the  initial  value  of  the  noise  filter  <^N  and  N  are  taken  as  oG, 
the  rma  of  noise.  °  0 

For  the  minidrone,  a  special  procedure  is  followed.  The  drone  is  assumed 
to  be  flying  at  some  load  factor  N^,  then  the  initial  angle-of-attack  and 
elevation  angle  are  computed  using  the  vehicle  aerodynamic  properties  as 
follows: 

Compute  the  dynamic  pressure 


using  the  initial  velocity  and  the  formulas  given  in  Section  2.3.  The  initial 
angle-of-attack  and  control  surface  angle  then  found  as 
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3.2  Update  all  Variables 

Updating  of  all  variables  consists  in  replacing  the  value  of  the  variable 
in  n+1  with  the  value  of  the  variable  in  n,  thus  for  a  typical  variable  x 


A  .  "  A 

n-1  n 

3.3  Guidance  Equations 

Using  the  block  diagrams  given  in  Section  2.5,  the  guidance  equations 
were  developed,  and  are  given  in  Tables  3.1  to  3. A.  In  considering  these 
tables,  it  is  noted  that  the  elevation  and  azimuth  gains  are  adjusted  in 
accordance  with  the  procedure  outlined  in  Section  2.5.  The  guidance  commands 
are  computed  in  the  sensor/processor  section  of  the 

Znrc  YMrF 

simulation  described  in  Section  3.7. 


3.4  Aerodynamic  Forces 

The  dynamic  pressure  and  Mach  no.  are  computed  as: 


p  p  (1  -  1.111  Z  +  0.36  Z  2) 
no  n  z 

p  -  0.002378 

o 

Zn  -  Zy  /40000 

n 

a  -  1117  -  (4.22  x  10-3)  Zy 

°n  n 

a  ■  965 

o 

n 


M 


n 


q 


n 


V  /a 
n  o 


0  <  Zy  <  36000 
n 

36000  <  Zy 

n 


Then  with  the  most  current  values  of  a,0,6_,6_,M,  the  aerody- 

n  n  e  k  n 

namlc  coefficients  are  computed  or  looked  up  in  a  table  depending  upon  the 
form  in  which  the  coefficients  are  provided.  The  form  for  each  of  the  vehicles 
is  given  in  Appendices  1-4. 
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TABLE  3.1 


MISSILE  GUIDANCE  EQUATIONS 
(PITCH) 


l*  AZC  “  ~°Z  XZMFF 
n  n 


"  AZBIAS  “  KG  (0GMI  +  eREF^* 

2.  LIMIT  A^c  to  +  20  g's  -  A£l(, 

n 

3.  LIMIT  Ag  /32.2  to  f  20  g'a  -  A^ 

n  n 

4*  *ZC  “  ^LC  “  ^L 


5*  xn  *  F2  xn  +  ^  F2^  KA  A ZC 


6. 


<l-^>  \  +  r  ka  *zc 

/  L  T1 


7.  LIMIT  <$'  to  +  14.14*  -  6' 


EL 


8.  LIMIT  Qn  to  100*/SEC  - 


9'  6E  m  F3  6E  .  +  (1  ~  V  (<SEL  +  KQ  QL  ^ 
n  n-x  n  n 


10.  LIMIT  6_  to  +  42.42*  -  6 


EL 


-T/t, 


-T/t, 


F2  -  e 


F3  -  e 


*Note  chat  this  term  Included  only  after 


9GM2 


<  9 

c 
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TABLE  3.1  (cont'd) 


1. 


MISSILE  GUIDANCE  EQUATIONS 
(YAW) 


“  -gy  \mff 

n 


2.  LIMIT  A^c  to  +  20  g's  -  A^CL 

n 

3.  LIMIT  Ay  /32 .3  to  +  20  g's  - 

n  n 

4.  Ayc  -  +  A^ 

n  n  n 

5‘  xn  ‘  F2  *„  +  <l  -  F2>  KA  *YC 

n 

6.  «  -  a  -  xn  +  £  KA 


7.  LIMIT  <SR  to  +  14.14" 
n 


6' 

RL 


8.  LIMIT  Rn  to  100° /SEC  -  R^ 


>•  '  F3  «R  ,  *  (1  -  F3>  <5Jl  +  h  h  > 

n  n-1  n  n 


10.  LIMIT  6r  to  +  42.42°  -  6^ 

n  ri 
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TABLE  3.2 


PROJECTILE  GUIDANCE  EQUATIONS 

A.  PITCH 

l'  *2Cn  *  ^ZMFP  -  *2.1*8  -  «G  (6OMl  +  W* 
2*  5e  “  ka  *zc 

n  n 

3.  LIMIT  61  to  +  12“  -  6' 

E  —  EL 

n  n 

4-  6En  -  SELn  +  Kq  «. 

5.  LIMIT  to  +  14°  -  6__ 

E  —  EL 

n  n 

B.  YAW 

1*  ^YC  "  ~Gy  ^YMFF 
n  n 


2. 

n 

ka\c 

n 

3. 

LIMIT 

\  “  ± 12  •  -  \ 

4. 

V 

5. 

LIMIT 

h  t0  ±  14’  ' 

*Note  that  this  term  is  included  only  when 


0 

GM2 


<  0 

c 


28 


TABLE  3.3 


MINIDRONE  GUIDANCE  EQUATIONS 


PITCH 

l'  4ECn  “  “GZtl  ^ZMFE  "  6EBUS  "  K0  <9GM1  +  eR£F)<+  KQ 

2'  4E  '  P3  V  ,  +  ‘  V  SEC 

n  n-1  n 

3.  LIMIT  6  to  +  15°  - 

E  EL 

n  n 

YAW 

X*  5RC  “  ~gy  Km?  +  ‘Si  Rn 

n  n 

2*  6R  "  F3  6R  .  +  (1  “  *3)  «RC 
n  n-1  n 

3.  LIMIT  <SR  to  +  15°  -  6^ 


*Note  that  this  term  13  included  only 

uh“  9  cm  9  ec 
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TABLE  3.4 

GLIDE  BOMB  GUIDANCE  EQUATIONS* 
(PITCH) 


ANGLE-OF-ATTACK  HOLD  (a-HOLD  SI,  S2,  S4-1) 
1.  CHECK  a 

IF  an  S.  «L1  xn  Mo‘“n  '  \i> 
IFan>L2  xn-Mo(an  -  a^) 

OTHERWISE  x  -  0 
n 

2-  5S  -  * 1 K,  +  Vi> 

n  n— l 


3.  61  -  6"  +  x 

a  £  XI 

n  n 


4.  LIMIT  6^  to  +  12°  -  5' 

n  n 

5-  6E  ’  6EL  +  ^  “  “REP  +  “n 


6.  LIMIT  6  to  +  15*  -  A 

E  EL 

n  n 


GIMBAL  HOLD  (S2,  S3-1) 


1.  CHECK  a 


IP  “»  -S  V  *„*  V°»  -  \l> 

IF  an  —  °L2  xn  "  Mo  <CV  -  “u> 

OTHERWISE  x  -  0 
n 


2'  fiE  "  Xn  "  KZH  (°GM1  +  6REp 
n 


3.  LIMIT  to  +  12°  -  6el 
n  ti 


*See  mode  sequencing  given  in  Table  2.3 
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TABLE  3.4  (cont'd) 


GLIDE  BOMB  GUIDANCE  EQUATIONS 
(YAW) 


A. 


<J>  COMMAND  (S7 ,  S8-1) 

3-  \  ■  Sh  -  V 

2.  LIMIT  5"  to  +  3°  ■  <S' 
R  -  R 


3. 

4. 


6  =6' 

R  RL  +  K_  »  R 
n  n  yB  n 

LIMIT  6  to  +  15°  -  6nT 
K  —  RL 


B.  ACCELERATION  AUTOPILOT  (S6,  S7-1) 

1*  Vf 

n  n 

2.  LIMIT  Ayc  to  NyL  • 

n  n 

3‘  SR  ■'yV  *  *Z  /32'2)  +  Kr  \ 


*•  \  (Sr'  +  SR  > 

o  n-I  n .  n-1 

5.  61  -  6"  + 


R„  +  iSfL  NcL 


n 


6.  LIMIT  61  to  +  38  «  6' 

R  -  RL 

n  n 


SR -sk  +Sn\ 

n  n 


8.  LIMIT  6r  to  +  15° 
n 


RRL 
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TABLE  3.4  (cont'd) 

GLIDE  BOMB  GUIDANCE  EQUATIONS 


*■  SE  + 


(YAW) 


LIMIT  6e  to  +  15°  -  6el 
n  n 


ACCELERATION  AUTOPILOT  (S2,  S4,  S5-1) 
1.  CHECK  a 


2. 

3. 

4. 

5. 

6. 

7. 

8. 


IF  “n  i  “LI  xn 

^  an  >  Ol2  *n 

OTHERWISE  x  -  0 
n 

AZC  "  “GZ  ^ZMFF 
n  n 

LIMIT  A^c  toNZL-AZCL 
n  n 

A’zc  “  KZ  (AZCL  ~  Az  /32.3  -  +  xq)  +  K 

n  n  n  ^ 


*'i  - 

s'i  *i 

+  AZC 

n 

n-1 

n 

n-* 

6'  - 

6"  +  x 

+  K„T 

A__ 

E 

n 

E  n 

n 

ZL 

ZCTj 

n 

LIMIT 

5E  * 
n 

5g  to  +  12°  - 
n 

fiEL  +  KZB<>n 
n 

6EL 

n 

LIMIT 

6*  to  + 
En 

15°  - 

Si 

-  Van-“L1> 

"  Mo<an  "  aL2) 


9. 


When  the  coefficients  have  been  obtained,  the  aerodynamic  forces  and 
moments  are  computed  as  follows: 

,L  ,  "»  S  ^  ^ 

\  \  Vv 

q  S?r 

M  -  C  -- 

Zn  n  lZZ 

n 

q  S 

A  -  C  — 

x  Y  m 
n  n  n 


\ 


q  S 

n  n 

“Y  m 
n  n 


q  S 

AZ  “  CZ  ~ 
n  n  n 

where  the  subscript  has  been  included  on  Iyy*  Igg  to  account  for  the  variation 
In  the  Inertial  parameters  caused  during  the  burn  period  of  powered  vehicles. 

3.5  Equations  of  Motion 

•  •  •  •  • 

The  angular  and  translational  accelerations  Q_ ,  R  ,  u  ,  v  ,  w  are  comp- 

n  n  n  n  n 

uted  from  the  equations  given  in  Section  2.2  using  the  most  current  values  of 
the  variables.  These  equations  are  then  integrated  as 


Q  JQ  ,  +  4  (Q  +  Q  ,) 

yn  vn-l  2  xn  xn-l 


\  '  Vi  + 1  «»  +  W 


R  -  R  .  +  4  <R  +  R  .  > 
n  n-1  2  n  n-1 

.  .  ^  .  T  c_n _  n-1 

vn  ’n-1  2  ^cos6  cos0  ,  ■ 

n  n-1 


u  ■  u  ,  +  -r-  (u  +  u  ,) 
n  n-1  2  n  n-1 
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V  -  V  .  +  ■=■  (v  +  V  ,) 
n  n-1  2  '  n  n-1 


W  -  w  ,+■=■(»  +  w  ,) 

n.  n-1  2  n  n-1 


3.6  Kinematics 

Having  the  velocity  components  in  body  axes  (un,  vQ,  wn,)  as  well  as  the 

orientation  of  the  body  axis  (0n,  fn);  the  velocity  of  the  vehicle  in  space 

X  ,  Y  ,  Z  ,  can  be  computed  from  the  transformation  equations  given  in 
vn  vn  n 

Section  2.1.  Finally  the  position  of  the  vehicle  in  space  is  obtained  by 
integration  as 


t  +  4  (X  +  X  ) 
v.  2  v  v. 
m-1  n  n-1 


and  the  same  for  the  Y  and  Zw  components. 

n  n 


3.7  LOS  Processor 

As  noted  in  Section  2.4,  the  measured  LOS  angles  are  given  by 


XZM  "  XZ  +  0eN  +  9eB  +  0eR 
n  n  n 


YM 


XY  +9CK  +9EB  +  *eR 


n  n  n 

where  Xz  and  XY  are  computed  from  the  equations  given  in  Section  2.4. 
n  n 

The  differential  equations  for  the  noise  components  of  the  body  rate 
measurements  must  be  cast  into  a  difference  form  for  use  in  the  digital  simul¬ 
ation.  The  white  noise  input  into  these  equations  is  approximated  by  the 
sequence  from  a  random  number  generator  as 


n  -  K(R  -  0.5) 
n  n 

where  is  a  random  number  between  0  and  1  and  K  is  a  constant  chosen  to  make 
the  spectral  density  of  the  sequence  unity.  For  a  sample  time  of  T,  the  spec- 
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trum  of  Che  sequence  nn  is  a  flat  double-sided  function  between  ±  idq  (the 
folding  frequency  given  as  (2ir/T)/2rad/sec) .  The  intensity  of  the  noise  is 
then 

c  _  K2  n  -  .  ir 

o  24oi  T  o  T 

o 

For  unit  intensity  (S0  *  1) 

Then,  in  difference  form 


QeN  “  A  <*eN  .  +  (1_A)  °G  2*G  nn 
n  n-1 

Rc»  -ARC»  ,  +  <1-*)  °G  “g  K 
n  n-1 


-T/tg  *  t 

where  A  »  e  and  nn  distinguishes  the  sequence  of  n  from  nn-  The  initial 

condition  of  these  filters  is  *  ReJJ  -  Op.  Then, 

o  o 


-e«  ,  +  T  +  <U  > 

n  n-1  o  n-1 


♦eB  *+W  .  +  T(ReN  /cot,6n  +  ReN  /c0«en-i) 

n  n-1  n  n-1 


For  the  gyro  bias  component 


1 

eB 

-  6 

eB  , 

+  T^eB 

n 

n-1 

0 

eB 

“^B  , 

T 

+  T  <\b 

n 

n-1 

and 

*>•„  "e 

the  gyro 

ation  run.  Finally,  for  wide  band  receiver  noise 
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%  -  m  \  <*.  -  °-5> 

n 

♦a  ‘  /ir \ <K  '  °-5> 

n 

i 

and  oR  is  Che  rms  value  of  the  receiver  noise  and,  again,  the  Rn  distinguishes 
the  random  sequence  for  <|»eR^  from  that  of  6gR^. 

The  next  step  is  to  filter,  differentiate  and  filter  as  noted  in  Section 
2.4  the  difference  equations  for  elevation 


"  *Z  XZMF  ,  +  (1_AZ*  (XZM  * 


XZMF  *  (XZMF  “  XZMF  )/At 
n  n  n-1 


Wf  *  ^  XZMFF  .  +  (XZMF  5 

n  n-1  n 


SjMFFF  "  h  XZMFFF  . 
n  n-1 


+  (1-^)  x2 


-Ltli . 


where  A  -  e 
z 


and  At  ■  t  -  t  , .  Similar  equations  are  used  for 
n  n-1 


azimuth.  The  values  of  the  LOS  rates  in  elevation  and  azimuth  are  used  to 
povlde  commands  to  the  guidance  computer. 

In  a  like  manner,  the  filtered  glmbal  angles  are  computed  as 


9GM1  “  AA  flGMl  +  (1_V  (XZM  "  V 
n  n-1  n 


0GM2  "  ^  8GM2  ,  +  (1_AB)  eGMl 


where  A^  ■  e  TA  and  A^  »  e  TB 
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4.0  SAMPLE  NUMERICAL  RESULTS 

Figures  4.1  to  4.4  present  sample  simulation  output  for  the  four  vehicle 
types  described  in  the  previous  sections.  These  runs  include  only  the  effects 
of  gyro  noise  and  bias.  Except  for  the  projectile  for  which  the  gyro  bias  was 
taken  as  0.1°/sec,  the  bias  value  was  0.06°/sec  and  the  noise  rms  was  10°/hr 
with  a  correlation  time  of  15  minutes  [12].  Only  the  bias  contributes 
significantly  to  the  error. 

As  may  be  seen,  Fig.  4.1(a)  shows  the  overall  trajectory  in  the  pitch  of 
elevation  plane.  The  gyro  bias  caused  the  vehicle  to  climb  shortly  after 
launch.  The  numbers  printed  on  the  trajectory  indicate  flight  time  after 
launch  in  seconds  (or  for  very  long  flights  in  the  case  of  the  minidrone  in 
Fig.  4.3(a),  in  minutes). 

The  projectile  trajectory  in  Fig.  4.2(a)  has  its  starting  condition  at  a 
height  and  speed  consistent  with  roughly  the  midpoint  of  its  flight  [9].  Since 
the  minidrone  is  a  loiter  vehicle,  its  starting  condition  was  approximately 
level  (1  g)  flight. 

The  starting  points  of  the  missile  and  glide  bomb  trajectories  (4.1(a) 
and  4.4(a),  respectively)  were  based  upon  typical  conditions  at  the  end  of 
midcourse. 

Figures  4.1(b),  4.2(b),  4.3(b),  and  4.4(b)  give  the  time  history  of  the 
vehicle  velocity,  the  elevation  acceleration  command  and  the  pitch  load  factor 
actually  achieved  by  the  vehicle.  In  the  case  of  the  missile  (fig.  4.1(b)), 
the  velocity  is  shown  increasing  from  about  800  ft/sec  at  the  start  to  a 
maximum  of  about  1360  ft/sec  after  4  seconds,  which  results  from  the  rocket 
propulsion  after  launch.  After  bum  out,  the  velocity  markedly  decreases  with 
time  until  a  fairly  stable  level  of  574  ft/sec  is  achieved.  The  initial  value 
of  the  commanded  load  factor  represents  the  bias  command  used  to  shape  the 
trajectory.  The  dip  in  the  pitch  load  factor  during  the  time  after  burn  out  is 
due  to  the  correcting  command  needed  to  bring  the  vehicle  back  on  course  from 
the  climbing  maneuver  caused  by  the  gyro  bias  error. 
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GROUND  RANGE  (nmi) 


Fig.  4.1(a).  Sample  output  for  missile 


PITCH  LOAD  VELOCITY 


PITCH  BOOT  ANGLE  PITCH  BOOT  RATE  ANGLE 

(deg)  (deg)  (deg) 


Fig.  4.1(c).  Sample  output  for  missile. 


40 


YAW  LOAD  FACTOR  YAW  LOAD  FACTOR  CROSSRANGE 

(g-units)  COMMAND  (g-units)  (left) 


Fig.  4.1(d).  Sample  output  for  missile. 
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AZIMUTH  ANGLE  YAW  BODY  RATE  ANGLE 

(deg)  (deg/s)  (deg) 


140  130  120  110  100  90  80  TO  60  50  40  30  20  10  0 

TIME  (s) 


Fig.  4.1(e).  Sample  output  for  missile 


ALTITUDE  (ft) 


4.500 
4,000 
3,000 

3.500 

2.500 

2,000 

1.500 

1,000 

500 

0 


PITCH  LOAD 

PITCH  LOAD  FACTOR  COMMAND  VELOCITY 

FACTOR  Ig-uniU) 


PITCH  BODY  ANGLE  PITCH  BODY  RATE 
(deg)  (deg/s) 


-10 
-20 
-30 
-40 

28  26  24  22  20  1 8  16  14  12  10  8  6  4  2  0 

TIME  (s) 

Fig.  4.2(c).  Sample  output  for  projectile. 
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YAW  LOAD  FACTOR  YAW  LOAD  FACTOR  CROSSRANGE 

(g-units)  COMMAND  (g-units)  (kft) 


AZIMUTH  BOOT  ANGLE  YAW  BODY  RATE  ANG 

(deg)  (deg/s)  (de 


Fig.  4.2(e).  Sample  output  for  projectile. 


Fig.  4.3(a).  Sample  output  for  minidrone. 


8 


PITCH  LOAD 

PITCH  LOAD  FACTOR  COMMAND  VELOCITY 

FACTOR  (g-unit*)  (g-units)  (ft/#) 


PITCH 


TIME  (s) 


Fig.  4.3(c).  Sample  output  for  minidrone. 


YAW  LOAD  FACTOR  YAW  LOAD  FACTOR  CROSSRANGE 

(g-unitt)  COMMAND  (g-unit*)  (kft) 


SIDESLIP  ANGLE 

YAW  ACHIEVED  CONTROL  SURF  DEFL 


04310-H 


GROUND  RANGE  (nmi) 


Fig.  4.4(a).  Sample  output  for  glide  bomb 


PITCH  BOOY 
ANGLE  (deg) 


-20  -  J 

-30  - 

-4°brTi _ i _ i _ i _ i  .  j _ i _ l _ i _ i  . .  i  ~ 

120  HO  100  90  80  70  60  50  40  30  20  10  0 

TIME  (•) 

Fig.  4.4(c).  Sample  output  for  glide  bomb. 
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YAW  LOAD  FACTOR  YAW  LOAD  FACTOR  CROSSRANGE 

(g-units)  COMMAND  (g-unit«)  (kft) 


Fig.  4.4(d).  Sample  output  for  glode  bomb. 
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AZIMUTH  BODY  YAW  BODY  ANGLE 

ANGLE  Ideg)  RATE  (deg/s)  (deg) 


TIME  (8) 

Fig.  4.4(e).  Sample  output  for  glide  bomb 


For  the  projectile  (4.2(b)),  there  is  much  less  variation  in  the  velocity 
from  the  starting  value.  The  decrease  indicates  a  drag  of  between  V4  and  V2  g 
which  is,  in  part,  balanced  by  the  gravitational  component.  After  about  20 
seconds,  the  gravitational  component  becomes  slightly  larger  than  the  drag 
component  and  the  vehicle  velocity  increases  slightly.  The  abrupt  change  in 
the  load  factor  command  at  the  beginning  of  the  flight  is  the  result  of  the 
vehicle  coming  to  equilibrium  when  its  control  surface  is  released  from  a  2ero 
fixed  position  at  the  start  of  the  flight.  The  command  is  a  combination  of  the 
several  Inputs  described  in  sections  2.5  and  3.3. 

For  the  minidrone  (Fig.  4.3(b)),  the  flight  is  fairly  uneventful  since  it 
was  assumed  to  be  in  straight  and  level  flight  at  the  start.  The  slight 
sagging  of  the  trajectory  is  caused  by  the  gyro  bias  error. 

For  the  glide  bomb  (Fig.  4.4(b)),  the  velocity  decays  gradually  from  the 
initial  condition  of  700  ft/sec  to  about  570  ft/sec  over  a  period  of  about  82 
seconds.  Then,  the  command  of  the  vehicle  is  switched  from  angle-of-attack 
hold  in  pitch  to  glmbal  angle  hold,  then  about  15  seconds  after  this  to  the 
terminal  mode  of  proportional  navigation  and  the  vehicle  begins  to  dive, 
picking  up  speed  until  the  final  terminal  velocity  of  about  900  ft/sec  is 
reached.  There  is  a  substantial  transient  associated  with  the  switching  of 
the  guidance  mode  as  evidenced  by  the  load  factor  response  also  shown  in  Fig. 
4.4(b).  It  is  appropriate  to  note  that,  although  it  is  included  on  the  center 
plot  of  4.4(b),  the  commanded  pitch  load  factor  is  not  used  until  the  terminal 
mode  is  switched  in  at  about  97  seconds.  As  may  be  seen,  the  command  level  is 
about  0.65g  at  this  time. 

Figures  4.1(c)  to  4.4(c)  provide  angle-of-attack  and  control  surface 
deflection  along  with  body  pitch  rate  and  angle.  The  transients  in  the  projec¬ 
tile  noted  above  are  particularly  noticeable  in  the  body  rate  of  Fig. 
4.2(c).  The  switching  transients  for  the  glide  bomb  vehicle  are  quite  notice¬ 
able  in  the  response  of  all  the  variables  in  Fig.  4.4(c). 

Figures  4.1(d),  4.2(d),  4.3(d),  and  4.4(d)  provide  the  cross  range  time 
history,  as  well  as  the  command  and  achieved  yaw  factor.  Figures  4.1(e), 
4.2(e),  4.3(e),  and  4.4(e)  give  the  slide-slip  angle  and  rudder  control  angle 
along  with  the  yaw  body  rate  and  azimuth  angle.  The  scheduled  reduction  in 
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Che  yaw  gain  during  Che  picch  gimbal  hold  and  terminal  mode  switching,  as 
discussed  in  Section  2.5,  can  be  clearly  seen  on  the  plots  of  Figures  4.4(d) 
and  (e) .  Mote  the  slight  transients  in  the  yaw  body  rate  due  to  the  switching 
as  well  as  the  abrupt  change  in  the  level  of  the  commanded  yaw  load  factor. 
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APPENDIX  1 


MISSILE  AERODYNAMIC  DATA 


The  data  presented  In  this  Appendix  is  based  upon  reference  5.  In  terms 
of  the  tabular  aerodynamic  data,  we  define  the  following  parameters: 


_N 

a 

s 

6 

C 

m 

a 


— !i  - 

B 

AC 


E 

AC 
_ n 

TT 


CNl  +  CN2 |aj  +  CN3  a2 


CN1  +  CN2  j  B |  +  CN3  B* 


CM1  +  CM2  |a|  +  CM3  a 


CM1  +  CM2 | $ |  +  CM3  82 


DCM1  +  DCM2  a 


DC Ml  +  DCM2  0 


and 


C  -  CXI  +  CX2(a2  +  82)  +  CX3(a2  +  g2)1/2  (6_2  +  <S_2)1/2  +  CX4(6  2  +  6  2) 


R 


'r "  ^  *(r) s 
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where 


C 

*Q 

Al 

AZ 

AZ 


•  C  -  -573 

*R 

*  0.1166  -  0.1964t 

-  0.0184  -  0.0362  (t  -  0.5) 

-  -0.1084 

-  3.5 

-  Cz  -  1.0 


0  <  t  <  0.5 
0.5  <  t4 
4  <  t 


and  Che  needed  aerodynamic  parameters  are  given  In  Table  A.l  as  functions  of 
Mach  No.  Note  that  a,  8,  <5  ,  and  6_  are  in  degrees. 

Since  the  Inertial  parameters  are  changing  during  the  first  four  seconds 
of  flight . (after  launch)  the  following  are  also  needed  for  the  equations  of 
motion: 


m  -  14.3478  -  1.3756t 

0  < 

t  <  0.5 

m  -  13.66  -  0.3714  (t  -  0.5) 

0.5 

<  t  <  4 

m  *  12.36 

A  < 

t 

IOT  -  -  64  -  3.8t 

0  < 

t  <  0.5 

YY  ZZ 

1^  -  Izz  -  62.1  -  1.0286  (t  -  0.5) 

0.5 

<  t  <  4 

hr  ’  hi  ’  58'5 

4  < 

t 

Tx/«  '  880 

0  < 

t  <  0.5 

Tx/»  '  160 

0.5 

<  t  £  4 

O 

1 

f3 

. 

•r 

4  < 

y 
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TABLE  A.l 

MISSILE  AERODYNAMIC  DATA 


M 

CXI 

CX2 

CX3 

CX4 

1 

0-00 

-3. 25000E--01 

1.80000E-04 

-1 . 20000E-03 

-  3.50000E-  04 

2 

0.30 

-  3. 25000E-01 

1 . 30000E  04 

-1.20000E-03 

-3.50000E-04 

3 

0.£0 

-3.25440E-01 

1.25540E -04 

- 1 . 08000E-03 

-3. 39000E04 

4 

0.80 

•3.27557E-01 

3.G91G0E-05 

-8 . 37000E -04 

-3. 10000E- 04 

5 

0.95 

■ 5.50570E-01 

1.21040E  04 

-9. 54000E-04 

-3. 11700E-04 

6 

1.  14 

-8.  73650E-01 

-1.21470E04 

-G. 31000E-04 

-  2. 00700E-04 

7 

1.30 

--3.48300E-01 

-2.04900E -04 

-G.G5500E-04 

-3. 72000E-  04 

8 

1.80 

- 1 . 02000E  +  00 

-2.60000E-04 

-&. 30000E-04 

-4. OOOOOE-  04 

3 

2.00 

- 1 . 09000E+00 

-3.20000E  04 

- 7 . 20000E -04 

- 5. OOOOOE- 04 

10 

2.50 

-1. OOOOOE +00 

-3.  20000E--04 

-7.20000E-04 

-5. OOOOOE  04 

M 

CN1 

CN2 

CN3 

1 

0.00 

-1.90000E-01 

-5. 08500E  03 

8. 25000E-05 

2 

0.30 

-1.90000E  01 

-5.08550E-03 

8. 25000E-05 

3 

O.GO 

-  1.32420E-01 

-G.G0250E-03 

1. 11400E-04 

4 

0.80 

- 1.34710E-01 

-S.63130E  03 

1.20690E-04 

5 

0.95 

-2.01370E-01 

-6. 59920E -04 

9. 08910E-05 

6 

1.14 

-2.073G0E01 

-£. 6G600E- 03 

1 . 01 500E-04 

7 

1.30 

-2. 10810E-01 

-G.341G0E03 

8. 52700E-05 

8 

1.80 

- 1.30000E01 

-5. 59500E  03 

9. 50000E-05 

9 

2.00 

•1 . 83000E-  01 

-5. 78000E-03 

1.01000E-04 

10 

2.50 

-  1.83000E-01 

-S.78000E  03 

1 . 01000E-04 

M 

CM1 

CM2 

CM3 

1 

0.00 

-3. 30000E- 02 

-8. 30000E  03 

1 . 25000E -04 

2 

0.30 

-3. 30000E-02 

-8. 30000E  03 

1.25000E  04 

3 

O.&O 

- 5. 33280E-02 

-B.1S32OE03 

2. 04750E-04 

4 

0.80 

~G. 7G040E-02 

-3.028G0E-03 

2. 51810E-04 

5 

0.35 

-7. 31390E-02 

-1. 13589E  02 

3. 87040E-04 

£ 

1.14 

-7.83710E-02 

-G. 51000E  03 

2. 47470E-04 

7 

1.30 

-7.00940E-02 

-3. 55190E- 03 

3. 1 3750E-04 

8 

1.80 

-4.47500E-02 

-3.22550E  03 

2. 70750E-04 

9 

2.00 

-2. OOOOOE-02 

-9. 570O0E- 03 

2.50000E-04 

10 

2.50 

-2. OOOOOE -02 

-3. 57000E  03 

2. 50000E-04 

M 

DCM1 

DC  M2 

1 

0.00 

- 1 . G2000E-01 

- 4 . 86000E  ~05 

2 

0.30 

-1 . G2000E-01 

-4.8GOOOE  05 

3 

0. 60 

-2. 10180E-01 

-2.  52900E-  05 

4 

0.80 

-  2. 32280E-01 

-8. 49B00E  OG 

5 

0.95 

-2. 125G0E-01 

-5.701G0E  05 

6 

1.  14 

-1.G4490E-01 

-2. 75700E*  05 

7 

1.30 

-1 . 44850E-01 

-2.&1600E  05 

8 

1.80 

-1 . 20000E-01 

-2.50000E  05 

9 

2.00 

-1 . 00000E-01 

-2.40000E  05 

10 

2.50 

-l.OOOOOE-Ol 

-2.  40000E-  05 

• 
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APPENDIX  2 


PROJECTILE  AERODYNAMIC  DATA 

The  data  presented  In  this  appendix  Is  based  on  references  6,  7,  and  8. 
In  terms  of  tabular  aerodynamic  data,  we  define  the  following  parameters: 


—  -  CNALPHA 
a 


CNALPHA 


C- 

~  -  CNDELTA  (M,  loti) 
6E 

°Y 

-A  -  -CNDELTA  (M,  |$|) 
°R 


C 

^  -  CMALPHA 


-CMALPHA 


-  CMDELTA  (M,  |a|) 
°E 


-A  -  CMDELTA  (M,  |g|) 
R 


C  .  -  -  [CXO  +  DELCX  (M,  a,  6_)  +  DELCX  (M,  B,  6_)] 
x  E  •  R 


c 

“r  ■  f  6  +  ^  5* 

c,  C 

c2  ■  r  a  +  i  ei 

C  ■  —  a  +  -A  S  +  ~  C  Q 

m  a  6_  E  r  2U  m.  y 

E  o  Q 

C  C  1  Cv 

Cn  “  7T  6  +  IT  6R  +  7“  Cn  R 

n  B  6R  R  rQ  2U  nR 


and  C L  ■  C  -  -1SS.  The  needed  aerodynamic  parameters  are  given  In  Table  A. 2 
Q  °R 

as  functions  of  M,  a,  B,  6  and  6  .  Note  that  DELCX  Is  a  three  dimensional 

E  R 

table.  All  of  the  angular  arguments  are  given  In  degrees. 
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TABLE  A. 2 


1 


j 


.  PROJECTILE  AERODYNAMIC  DATA 

M 

CXO 

CNALPHA 

CMALPHA 

1 

0.40 

0.355 

-0.  340 

-0.310 

2 

0.50 

0.  340 

~0. 343 

-0.313 

3 

0.60 

0.340 

-  0. 348 

-0.318 

t* 

0.70 

0.350 

-0.  360 

-0.323 

s 

0.80 

0.375 

-0. 385 

-0.335 

6 

0.90 

0.500 

-0.425 

-0.360 

7 

1.00 

0.670 

-0.440 

-0. 385 

CNOELTA 

FUNCTION 

OF  W  & 

ADS  VAL  AOA 

M 

A0A*0 

A0Aa5 

AOA®  1 0 

AOA= 1 5 

AOA«20 

1 

0.40 

-0. 108 

-0.110 

-0. 108 

-0. 045 

-0.023 

2 

0.50 

-0. 110 

-0.115 

-0.110 

-0. 048 

-0.025 

3 

0.60 

-0. 112 

-0. 120 

-0.112 

-0. 050 

-0.028 

4 

0.70 

-0.117 

-0. 125 

-0.117 

-0.055 

-0.030 

5 

0.80 

-0. 120 

-0. 130 

-0. 120 

-0.060 

-0.035 

6 

0.90 

-0.130 

-0. 140 

-0.130 

-0.075 

-0.050 

7 

1.00 

-O. 140 

-0.148 

-0.140 

-0. 100 

-0. 080 

CMDELTA 

FUNCTION 

OF  M  & 

ADS  VAL 

AOA 

M 

A0A»0 

A0Aa5 

AOA® 10 

AOA® 15 

AOA =20 

1 

0.40 

-0.340 

-0. 330 

-0.325 

-0. 120 

-0.110 

2 

0.50 

-  0. 348 

-0.  338 

-0. 328 

-0.110 

-0. 100 

3 

0.60 

-0.358 

-0.  345 

-0. 330 

-  0. 100 

-0.088 

4 

0.70 

-0.370 

-0. 357 

-0.340 

-0. 095 

-0.075 

5 

0.80 

-0. 386 

-0. 370 

-0. 350 

-  0. 105 

-0.073 

6 

0.90 

-0.413 

-0. 395 

-0. 365 

-0. 130 

-0.088 

7 

1.00 

-0. 445 

0.  422 

-0.383 

-  0. 163 

-0. 122 
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TABLE  A. 2  (cont’d) 
PROJECTILE  AERODYNAMIC  DATA 


DELCX  FOR  M=0.4 
CONTROL  SURFACE  DEFLECTION 


ADA 

-20 

-  15 

-10 

-5 

0 

5 

10 

15 

1 

20 

0.23 

0.28 

0.  IS 

0.  10 

0.02 

-  0.07 

-0.  IS 

-0.  19 

a 

-15 

0.28 

0.28 

0.  13 

0.  10 

0.02 

-0.07 

-0.09 

-0.03 

3 

-10 

0.28 

0.28 

0.  18 

0.  10 

0.02 

0.00 

0.03 

0.  13 

4 

-5 

0.24 

0.24 

0.  16 

0.03 

0.00 

0.02 

0.  10 

0.21 

5 

O 

0.52 

0.  30 

O.  15 

0.06 

0.00 

0.06 

0.  18 

0.22 

6 

5 

0.41 

0.21 

0.  10 

0.02 

0.00 

0.03 

0.  1G 

0.24 

7 

10 

0.26 

0.  13 

0.03 

0.00 

0.02 

0.  10 

0.  18 

0.28 

8 

15 

0.05 

-0.03 

-0.09 

-0.07 

0.02 

0.  10 

0.  18 

0.28 

9 

00 

-0.  10 

-0.  13 

-0.  18 

-0.07 

0.02 

0.  10 

0.  13 

0.  37 

DELCX  FOR  M=0. 8 
CONTROL  SURFACE  DEFLECTION 


AOA 

-20 

-  15 

-10 

-5 

0 

5 

10 

15 

1 

™20 

0.37 

0.25 

0.  10 

0.00 

-0.  10 

-0.  18 

-0.20 

0.37 

2 

-  15 

0.37 

0.25 

0.  10 

0.00 

-0.08 

-0.  1 1 

-0.05 

0.23 

3 

-10 

0.23 

0.20 

0.09 

0.00 

-0.02 

0.01 

0. 12 

0.23 

4 

-5 

0.29 

0.20 

0.09 

0.00 

0.03 

0.  11 

0.25 

0.40 

5 

0 

0.30 

0.  16 

0.07 

0.00 

0.05 

0.  16 

0.22 

0.37 

6 

5 

0.25 

0.11 

0.03 

0.00 

0.09 

0.20 

0.29 

0.30 

7 

10 

0.  12 

0.01 

-0.02 

0.00 

0.03 

0.20 

0.29 

0.09 

8 

15 

-0.05 

-0.11 

-  0.  08 

0.00 

0. 10 

0.25 

0.37 

-0. 15 

9 

20 

-0.20 

-0.  18 

-0.  10 

0.00 

0. 10 

0.25 

0.37 

0.45 

DELCX  FOR  M=1 . 0 
CONTROL  SURFACE  DEFLECTION 


AOA 

-20 

-15 

-10 

-5 

0 

5 

10 

15 

1 

—  20 

0.45 

0.35 

0.20 

0.03 

-0.  13 

-  0. 16 

-0. 13 

0.40 

2 

-15 

0.40 

0.29 

0.15 

0.02 

-0. 05 

-0.05 

0.03 

0.  35 

3 

-10 

0.35 

0.24 

0.08 

0.01 

0.01 

0.05 

0. 16 

0.30 

4 

-5 

0.30 

0.  13 

0.07 

0.00 

0.03 

0.  10 

0. 10 

0.25 

5 

0 

0.60 

0.40 

0.20 

0.03 

0.00 

0.06 

0. 18 

0.28 

6 

5 

0.42 

0.25 

0.  10 

0.03 

0.00 

0.07 

0. 13 

0.30 

7 

10 

0.30 

0. 16 

0.05 

0.01 

0.01 

0.08 

0.24 

0.35 

8 

15 

O 

■ 

4> 

0.03 

-0.05 

-0.05 

0.02 

0. 15 

0.23 

0.40 

9 

20 

-0.06 

-0.  13 

-0.  16 

-0.13 

0.03 

0.20 

0.  35 

0.45 

y 
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V 


°L 


where  all  angles  are  In  degrees.  The  coefficients  are  given  as  follows: 


C  -  -0.06 
m 


C  -  -0.0300 
m 


C  -  -0.0223 
m<5 


C  -  -17.0 
"Q 

CD  •  0.048 
o 


K  -  19.09 
a 

Ke  -  6.03 

CL  -  0.34 
o 

-  0.0960 

a 

C.  -  0.0062 
L6 


C  -  0.0140 


C  -  -0.0140 
n6 


C  -  -10.0 
n_ 


Cl  -  0.0262 

e 

-  -0.0044 

o 


The  variation  of  the  thrust  with  velocity  is  given  as  [8]: 


400 
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APPENDIX  4 


GLIDE  BOMB  AERODYNAMIC  DATA 


The  data  presented  In  this  appendix  la  based  upon  references  10 ,11,  12. 
The  aerodynamic  coefficients  are  expressed  as 

Cx  -  -{CX0(M,  a)  +  [CXBETA+2(M,0>  02  +  CXDELt2 (M)  (6£2  +  <SR2)]) 

o 

-  C,  0  +  CYDEL(M)6r 


Cz  “  cz  +  cz  a  +  czdel<m>5e 

o  a 


C  -  C  (M,a)  +  DELCM(M,a| B| )  +  CMDEL(M)6_ 


c  -  DELC  (M,a, |B|)  SGN(B)  +  C  DEL(H)5 
n  n  n  K 


and  Cz  ■  -0.25 
o 


-0.10 


Cy  “  -0.029 

0 

Table  A.3  gives  the  aerodynamic  parameters  in  these  equations  as  functions  of 
Mach  No.  and  angles  of  attack  and  sideslip  (a  and  0)  as  noted  in  the  equations. 
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TABLE  A. 3 


GLIDE  BOMB  AERODYNAMIC  DATA 


M 

CMDEL 

CZOEL 

CnDEL 

CYDEL 

1 

0.  SO 

-0. 043000 

-0.00700 0 

-0.006300 

0. 014S000 

2 

0.70 

-0.038400 

-0-007900 

-0.006100 

0.0148000 

3 

0.90 

-0.040100 

-  0. 009800 

-0. 006800 

0.0140000 

4 

1.  10 

-0.031400 

-0.008000 

-0.004400 

0.0105000 

5 

.1 . 30 

-0.034400 

-0.006.S00 

-0.003100 

0.0079000 

CttiO  AS  A 

FUNCTION  OF  M  &  ALPHA 

M 

ALP=- 15 

ALP*  -10 

ALP=-S 

ALP=0 

ALP  =  5 

ALP= 1 0 

1 

0.50 

0.003000 

0.068000 

0. 060 SOO 

0.0343000 

0. 009300 

0. 042300 

0.70 

0.070000 

0.080000 

0.058400 

0.0306000 

0.010100 

0.054500 

3 

0.  '30 

0. 170000 

0. 183000 

0. 102300 

0,0150000 

0.081800 

-  0. 079500 

4 

1.  10 

0. 371000 

0.090000 

0. 150500 

-0.0200000 

•0. 315000 

-  0. 566000 

5 

1. 30 

0.230000 

0.0S9000 

0.079200 

-0.0920000 

-0.320500 

-0. S83000 

CXO  AS  A 

FUNCTION  OF  M  &  ALPHA 

M 

ALP  ■-■15 

ALP =-10 

AlP=- 5 

ALP=0 

ALP  =  5 

ALP  =10 

1 

0.  SO 

0.033000 

0.003000 

0.039800 

0.0566000 

0.008200 

0. 06.2800 

2 

0.70 

0.044000 

0.034000 

0.040400 

0.0570000 

0.008200 

0.047800 

3 

0.90 

0. 069000 

0.048500 

0.060100 

0. 0855000 

0.086100 

0. 082000 

4 

1.10 

0. 153000 

0, 159800 

0.154700 

0. 1248000 

0. 142600 

0. 151500 

S 

1.30 

0. 1 20000 

0.135900 

0. 149100 

0. 1583000 

0. 162000 

O.  160100 

C.X 

BETA*  2 

AS  A  FUNCTION  OF  M 

&  ALPHA 

M 

ALP  =  15 

ALP  =-10 

ALP  =5 

ALP=0 

ALP =5 

ALP  =10 

1 

0.  SO 

0.050000 

0.050000 

0.154000 

-0. 7360000 

-0. 122000 

-  0. 120000 

2 

0.70 

0.070000 

0.070000 

-0. 120000 

-0.9100000 

-0. 168000 

-0. 163000 

3 

0,90 

0, 560000 

0. 560000 

-0.320000 

-0.6800000 

-0. 370000 

-0. 370000 

4 

1.  10 

-0,070000 

-0.070000 

-0.246000 

-0. 3800000 

-  0. 290000 

-  0. 290000 

5 

1.30 

0.000000 

0.000000 

-0.072000 

-0.2000000 

-0. 188000 

-  0. 138000 

CX 

DO.  *3 

M 

0.50 

AS  A  FUNCTION  OF  M 

1 

0.51S000 

2 

0.70 

0. 740000 

3 

0.90 

1.069000 

4 

1.10 

0.837000 

5 

1.30 

0.670000 

TABLE  A. 3  (cont’d) 

GLIDE  BOMB  AERODYNAMIC  DATA 


DELCM  os 

M 

&  AOS(BCTA)  TOR  ALPHA*- 15 

M 

BETA*0 

BETA=4 

BETA* 8 

1 

0.50 

0.000000 

0.019300 

-0.004200 

2 

0.70 

0.000000 

-0.012800 

-0.027500 

3 

0.90 

0.000000 

-0. 035800 

-0.210600 

4 

1.  10 

0.000000 

■  0.036900 

-0.090700 

5 

1.30 

0.000000 

-0.013100 

-0.040700 

DELCM  VS 

M 

&  ADS (BETA)  FOR  ALPHA*- 10 

M 

BETA=0 

BETA=4 

BETA=3 

1 

0.50 

0.000000 

0.023000 

0.039700 

2 

0.70 

0.000000 

0.008400 

0.026600 

3 

0.90 

0.000000 

0 . 004000 

-0.049200 

4 

1.10 

0.000000 

-0.007500 

-0.009900 

5 

1. 30 

0.000000 

-0.002900 

0.006700 

DELCM  US 

M 

&  ABS(BETA)  FOR  ALPHA*- 5 

M 

BETA=0 

BETA=4 

BETA=3 

1 

0.  50 

0.000000 

0.026700 

0.083600 

2 

0.70 

0. OOOOOO 

0.030000 

0.080600 

3 

0.90 

0.000000 

0.043800 

0.112200 

4 

1.10 

0.000000 

0.021900 

0. 070900 

5 

1.30 

0.000000 

0.007300 

0.054100 

DELCM  VS 

M 

&  ADS (BETA)  FOR  ALPHA*  0 

M 

BETA=0 

BETA=4 

BETA=8 

1 

o.  SO 

0. OOOOOO 

0.018800 

0. 076000 

2 

0.70 

0.000000 

0.023800 

0.088100 

3 

0.90 

0.000000 

0.021000 

0. 125800 

4 

1.10 

0.000000 

0.041700 

0. 131800 

5 

1.30 

0.000000 

0.007300 

0.017500 

DELCM  VS 

M 

&  ADS (BETA)  FOR  ALPHA*  5 

M 

BETA*0 

BETA*4 

BETA=3 

1 

0.50 

0.000000 

0.017100 

0.039900 

2 

0.70 

0,000000 

0.019000 

0.050100 

3 

0.90 

0.000000 

0.059500 

0. 129800 

4 

1.10 

0.000000 

0.031300 

0. OSS 300 

5 

1.30 

0.000000 

-0.004200 

-0.041000 

DELCM  VS 

M 

&  AOS (BETA)  FDR  ALPHA*  10 

M 

BETA*0 

BETA=4 

BETA*8 

1 

0.50 

0.000000 

0.010800 

-0.004300 

2 

0.70 

0.000000 

0.022500 

-0.005700 

3 

0.90 

0.000000 

0.077300 

0.127900 

4 

1.10 

0.000000 

0.026100 

0.064800 

5 

1.30 

0.000000 

0.036400 

0.111800 

OCT  A*  12 
0.0278000 
-0.0132000 
0. 1684000 
-  0.0926000 
0.0306000 


BCTA=  12 
0.0439000 
0.0181000 
-0.0696000 
0.01 63000 
0.0189000 


BETA* 12 

0.0600000 

0.0494000 

0.0494000 

0.0925000 

0.0684000 


BETA* 12 
O. 0960000 
O. 1027000 
O. 1427000 
O. 1 662000 
0.0463000 


BETA*  12 
0.0618000 
O. 1067000 
0.2102000 
0.  1964000 
-0.0034000 


BETA* 12 

0.0225000 

0.0699000 

0.2542000 

0.2339000 

O. 1583000 


TABLE  A. 3  (cont'd) 

GLIDE  BOMB  AERODYNAMIC  DATA 

DELCn  VS 

&  AOS (BETA)  FOR  ALPHA*  15 

M 

BET  A=0 

BETA=4 

BETA=3 

1 

0.  50 

0. 000000 

0. 003000 

0.001000 

2 

0,70 

0. 000(300 

0.003000 

0.002600 

3 

0.  90 

0.000000 

0.010000 

0,006.000 

4 

1 .  10 

0 . 000000 

0,005000 

0. 006000 

5 

1,30 

0.000000 

0. 002000 

0. 003000 

DEI. 

_Cn  VS 

m 

&  ADS (BETA)  FOR  ALPHA*- 10 

M 

BETA=0 

BETA=4 

BETA=S 

1 

0.  50 

0.000000 

0. 005000 

0.003500 

2 

0.  70 

0. 000000 

0. 005000 

0. 006700 

3 

0.90 

0.000000 

0.01 3000 

0,011000 

4 

1 ,  10 

0.000000 

0. 005500 

0.009000 

5 

1.30 

0. 000000 

0. 002500 

0. 002500 

DELCn  VS 

m 

h  ADS (BETA)  FOR  ALPHA* -5 

M 

BETA=0 

QETA=4 

BETA* 8 

1 

0.50 

0.000000 

0.006500 

0.003000 

2 

0.  70 

0.000000 

0. 007800 

0, 009600 

3 

0,90 

0 , 000000 

0.015500 

0.019000 

4 

1.  10 

0.000000 

0.00S500 

0.011000 

r, 

1.30 

0.000000 

0. 006000 

0,010000 

DELCn  VS 

M 

&  ADS (BETA)  FOR  ALPHA*  0 

M 

BETA=0 

BETA*4 

BETA=3 

1 

0.  50 

0.000000 

0.007000 

0.012500 

2 

0.70 

0.000000 

0,003000 

0.015000 

3 

0,90 

0.000000 

0.012000 

0.021000 

4 

1.  10 

0.000000 

0.015000 

0.015000 

5 

1.30 

0.000000 

0.005500 

0,014500 

DELCn  VS 

&  AOS (BETA)  FDR  ALPHA*  5 

M 

BET  A*0 

BET A* 4 

BETA=3 

1 

0.  50 

0.000000 

0.009500 

0.021000 

2 

0.70 

0.000000 

0.013000 

0,021000 

3 

0.90 

0.000000 

0.013500 

0,024000 

4 

1.  10 

0.000000 

0.013000 

0.024000 

5 

1.30 

0.000000 

0.004500 

0.014500 

DELCn  VS 

M 

&  ADS (BETA)  FOR  ALPHA*  10 

M 

BETA=0 

BETA=4 

BETA=3 

1 

0.30 

0. 000000 

0.010300 

0.024300 

2 

0.70 

0.000000 

0.010400 

0.026500 

3 

0.90 

0.000000 

0.013500 

0.024000 

4 

1.10 

0.000000 

0,016800 

0.030700 

5 

1.30 

0.000000 

0.004500 

0.014500 

BETA* 12 

O. 0000000 

0.0004000 

0.0030000 

0.0100000 

0.0000000 


BETA*  12 
0. 0000000 
0. 0060000 
0.0110000 
0.0150000 
0. 0005000 


BETA*  12 
0.0130000 
0.0062000 
0.0210000 
0.0160000 
0.0130000 


BETA*  12 
0.0205000 
O. 0330000 
0.0410000 
0.0250000 
O. 0275000 


BETA*  12 

0.0355000 

0.0301000 

0.0420000 

0.0340000 

0.0205000 


BETA=12 

0.0430000 

0 . 0428000 

0.0420000 

0.0423000 

0.0205000 
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